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ABSTRACT 



We present data from Integral Field Spectroscopy for 3 supercompact UV-Luminous 
Galaxies (ScUVLGs). As nearby (z ~ 0.2), compact (R50 ~l-2 kpc), bright Paschen-a 
sources, with unusually high star formation rates (SFR=3— 100 Mq yr~^), ScUVLGs 
are an ideal population for studying detailed kinematics and dynamics in actively star- 
forming galaxies. In addition, ScUVLGs appear to be excellent analogs to high redshift 
Lyman Break Galaxies (LBGs) and our results may offer additional insight into the 
dynamics of LBGs. Previous work by our team has shown that the morphologies of 
these galaxies exhibit tidal features and companions, and in this study we find that the 
dynamics of ScUVLGs are dominated by disturbed kinematics of the emission line gas- 
suggestive that these galaxies have undergone recent feedback, interactions or merg- 
ers. While 2 of the 3 galaxies do display rotation, v/cr <1 - suggesting dispersion 
dominated kinematics rather than smooth rotation. We also simulate how these ob- 
servations would appear at z ^2. Lower resolution and loss of low surface brightness 
features causes some apparent discrepancies between the low-z (observed) and high-z 
(simulated) interpretations and quantitatively gives different values for v/a, yet simu- 
lations of these low-z analogs manage to detect the brightest regions well and resemble 
actual high-z observations of LBGs. 



Introduction 



Across cosmic time, processes driving galaxy formation are e xpected to traiisitioii from hierar- 
chical merging at early times to secular evolution at low redshift (jWhite &: Reeall978l : lAbadi et al 
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20031 : iKormendv &: Fisheijl2005l ). In the case of hierarchical merging, galaxies grow as halos merge, 
causing disturbed kinematics and dynamics. Secular evolution describes slow and steady mecha- 
nisms, including gas accretion and internal processes, and is consistent with smooth rotation. 

The Galaxy Evolution Explorer {GALEX) has uncovered a population of local {z = 0.1 — 
0.3), intensely star- forming galaxies. These ultraviolet luminous galaxies (UVLGs) are selected 
based on their Far-UV (FUV) luminosities {J-ifuv > 2xlO^°L0). UVLGs exhibiting the highest 
surface-brightnesses (Jfuv > 10^ Lq kpc~^), named supercompact UVLGs (ScUVLGs), are ideal 
candidates for Integral-field Spectroscopy (IFS). They are bright sources whose compact sizes (with 
half-light radii ~ 1 — 2 kpc) accommodate the limited fields-of-view (FOV) of many IFS instruments, 
typical redshifts place the Paschen-a (Pa-a) line conveniently in the K-band filter, and unusually 
high star-formation rates (SFRs), with SFR=3— 100 Mq yr~^, designate them as a unique sample. 
These ScUVLGs bear striking resemblances to Lyman Bre ak Galaxies (LBGs) - s haring similar 
specif ic SFRs, metallicities, k inematics, and attenuations ( Heckman et al. (j2005l ). Hoopes et al 



(|2007l ). iBasu-Zych et al.l (120071^). Whe n resimulated at high redshift, their morphologies are similar 
to those of LBGs (lOverzier et al.ll2008l ). ScUVLGs constitute a local (z~0.2) sample well-suited for 
studying details of galaxy formation with very high physical resolution and sensitivity using IFS. 



LBGs, selected by the Lyman break dropout technique (ISteidel et al.ll200d'). f orm stars at 



Bouwens et al 



JSB) find only a 



intense rates, dominating the UV luminosity density at z >2. 
modest decrease in the UV luminosity density out to z =6, indicating that LBGs represent a major 
phase in early stages of galaxy formation and evolution. Many studies have used IFS to observe 
the dynamics of L BGs, at z = 2 — 3, in order to pinpoint the dominant formation pro c ess fo r 
these galaxies fi.e.. Law et all (l2007l. bood') . iNesvadba et al.1 (I2OO6I1 . IPorster Schreiber et al.l i^20od ). 



Wright et al.l tOOTi ) 



Nesvadba et al 



(|200a), and 



Stark et al 



((20081) )• These studies probe key 
questions: What causes such intense star formation in these systems - disk instabilities or tidal 
shocks from mergers? How does angular momentum evolve in early star-forming galaxies? Yet, at 
the high redshifts of LBGs, spatial resolution and depth are fundamentally limited. Except in the 
fortuitous but rare cases of gravitationally lensed LBGs, IF S at z >2 may be biased towards brighte r 
and larger LBGs whose dynamics differ from typical LBGs (jNesvadba et al.ll2008l : IStark et al.ll2008l ). 

We have observed 3 ScUVLGs with the OSIRIS instrumenllfl. While a larger sample of ScU- 
VLGs has been commissioned for IFS studies, in this Letter we share initial results and compare 
them with high- 2; LBG studies. In ^ we discuss data reduction, and analyze the results from these 
observations in ^ We discuss the relevance of this work, related to other high-z LBG studies in 
the final section, ^ Throughout our analysis, we apply ACDM cosmology with (Hq, ^m, O,a)={70 
km s-^Mpc-\ 0.3, 0.7). 



The data presented herein were obtained at the W.M. Keck Observatory, which is operated as a scientific part- 
nership among the Cahfornia Institute of Technology, the University of Cahfornia and the National Aeronautics and 
Space Administration. The Observatory was made possible by the generous financial support of the W.M. Keck 
Foundation. 
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2. Data and Analysis 



Using the OSIRIS (OH Suppressing InfraRed Imaging Spectrograph: iLarkin et al.l (j2006l )) in- 
strument at the Keckll telescope, we observed SDSS J092600. 40+442736.1 and SDSS J143417.15+020742.5 
(referred to as ScUVLG092600 and ScUVLG143417, respectively) in Spring, 2008 and SDSS J210358.74-072802.4 
(ScUVLG210358, henceforth) in Fa ll 2008 OSIRIS utihze s the L GSAO system (Laser Guide Star 
Adaptive Optics: IWizinowich et alJ (120061 ) : Ivan Dam et al.1 (120061 )) to improve spatial resolution. 



Several criteria were applied in selecting the targets and observing configuration. Using the 
LGSAO requires that targets be within 60" of a bright (M/j < 17 mag) tip-tilt reference star. 
Certai n filter and spaxal scale c ombinations result in excess thermal background and/or detector 
noise (jLaw. Steidel. &: Erbll2006l ). To optimize the sensitivity of our observations, we chose the Kn3 
filter with 50 mas lenslet"^ scale. For this combination, the resulting FOV for our observations was 
approximately 2" x 3" (the actual FOV varied because of slight re-centering during observations - 
see magenta boxes in Fig. [T]). The resolution and summary of our observations for each target are 
given in Tabled! Our velocity resolution is approximately 34-35 km s~^. Total on-target integration 
times were ~ 2 hours per target. 

As previously mentioned, selection of ScUVLGs depends on FUV luminosity and surface- 
brightness exceeding 10^'' '^ Lq and 10^ Lq kpc~^, res pectively. We exc l ude AGN from thi s samp le 
using two separate methods: r adio luminosit i es (se e 
and optical emission lines (see Hoopes et al. (j2007l ): 
ScUVLG sample consists of compact, starburst dominated galaxies exhibiting signs of small-scale 
merge rs or interactions on ly detected in the rest-frame optical at the resolution and sensitivity of 



Basu-Zvch et al.l moi): 



Best 



Baldwin. Philhps. fc TerlevichI (|l98lh ). The 



et al 



(12005!)) 



rpl 

i)) 



HST (lOverzier et al.l l2008l ) . Our 3 galaxies span stellar mass and SFR ranges of the full set of 



ScUVLGs. In Table (H we list some of their physical properties. Stellar masses are determined by 
fitting the spect ral energy d i stribu tion (SED) using unresolved Sloan Digit a.1 Sky Survey (SDSS) 
photometry (see lSalim et al.l (|2005l ) and lKauffmann et al.l (|2003l )'). assuming a lKroup al (j200ll ) Initial 
Mass Function (IMF). The SFRs are unattenuated SFRs, determined from radio luminosities (see 
Basu-Zvch et al.l (|2007l )). We estimate that typical uncertainties associated with the stellar mass 
and SFRs are ~30% based on comparing different mass modeling methods and SFR indicators. 
Average v alues and standard dev iations for these properties in the full sample of 45 ScUVLGs 
(studied in lBasu-Zvch et al.l hoQ^ )) are as follows: (SFR) = 21 Mq yr"^ and (Jsfr = 16 Mq yr'^ 
(log M^/Mq ) = 10.0 and cTiogM./Mo = 0.5. While both ScUVLG143417 and ScUVLG210358 are 
dominated by luminous, compact UV and line emitting regions, both contain significant low surface 
brightness disk structure and are more massive than the typical ScUVLG. 

We create 3-dimensional data cubes of spatial versus spectral dimensions by using the Keck/OSIRIS 
data reduction pipeline. The pipeline performs the following tasks: corrects detector non-linearity 
to form a data cube with calibrated wavelength units, subtracts a sky frame from the science frame 
to combine separate exposures into a single mosaicked data cube, removes cosmic rays, and cor- 
rects for detector cross talk. A custom IDL routine determines a mask for good data based on 
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a signal-to-noise (S/N) threshold set at Scr. After boxcar smoothing by a kernel of 3 lenslets in 
the spatial dimensions to increase the S/N, gaussians were fit to the spectral data to determine 
the peak, peak location, and FWHM of the Pa-a feature - creating S/N maps, velocity maps and 
velocity dispersion (a) maps, respectively. Resulting images for the 3 targets are shown in Fig. [2] 
(with Pa-a S /N contours marked for reference) . Fig. [1] displays the orientation and FOV of these 
targets. 

For ScUVLG092600, we also made po sition-velocity (PV) diagrams (shown in Fig. [3]) using 



the kpvslice tool in KARMA (|Goochlll996l l. An integrated spectrum combining all the unmasked 
data was created using custom IDL software. To calculate inclination, i, we assume that face-on 
(i =0) disks have circular isophotes and edge-on disks have intrinsic axis ratios of 0.20. Using Pa-a 
S/N maps, we use ratios of minor to major axes of the isophotes to measure inclinations. We use 
several isophote levels to determine the uncertainty of this measurement. 



Results 



HST images (shown in Fig. [T], see lOverzier et al.l ()2008l . l2009l )) reveal companions and disturbed 
morphologies, such as tidal features and plumes. Using IFS, we analyze the velocity structure for 
signs of smooth rotation. In two cases, we see evidence of regular rotation (ScUVLG092600 and 
ScUVLG210358), while ScUVLG143417 has no clearly defined velocity structure. However, even 
in the cases of ordered velocity, we find high dispersion (v/u <1), suggesting dispersion dominated 
systems rather than smoothly rotating thin disks. 



3.1. ScUVLG092600 

The bright component in the NW of this system has weak, but regular velocity structure 
across its disk. HST Ha imaging reveals a filamentary component extending into the south-west, 
potentially suggesting an outflow. Since this feature is undetected in optical or UV morphologies, 
it is unlikely to be a tidal feature or spiral arm. The mild velocity gradient may arise from an 
outflow in the south-western region of the disk, rather than from smooth rotation. For a rotating 
disk scenario, we find an inclination, i = 40 it 5°. The bottom panel of Fig. [3] displays the PV 
diagram, sliced along the central region of the disk (shown in magenta on the velocity field, at top). 
We determine a circular velocity of Vc = ^^^^i — ~3ll3 km s~^ . Qualitatively, the PV diagram 
shows that velocity dispersion dominates the velocity structure across the disk. Quantitatively, the 
integrated spectrum gives a =84 ±3 km s^^, and v/o" for this component is 0.37ib0.05. Assuming 
that the observed velocity dispersion of the ionized gas for this component can provide an estimate 



^While typical spiral galaxies have intrinsic axis ratios (qo) of 0.2, ScUVLGs may be thicker. We note that even 
for qo=0.5, values of Vc and v/a differ by < 10%, within measurement uncertainties. 
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for the disk's dynamical mass, M^y„, we use the following relation: 

Ca'^R 

Mdyn = (1) 

where C is the geometric constant (C = 3.4 for thin rotating disks, C = 5 for spherical geometry), G 
is the gravitational constant, R is the radius of the disk, and a is the velocity dispersion. Adopting 
the following values: R=0.3"=0.9 kpc, a =84 km s~\ we determine Mdyn =4.4-6.5x10^ Mq, 
corresponding to 3.4<C<5. 

The companion in the south-east, spanning few velocity resolution units across, has systemic 
velocity of = -|-80 km s~^ compared to the brighter component. While we do not detect any 
ionized gas linking these two components, they are likely in the process of interacting given their 



close proximity (< 3 kpc). As further evidence of interaction, iGoncalves et al.l (|2009l ) find that 
the bright component's central region is consistent with regular rotation, while outer edges become 
more disturbed and dispersion dominated. 



3.2. ScUVLG143417 

This galaxy's morphology differs from ScUVLG092600, with two bright cores (separated by 
<1 kpc) embedded in a fainter, disturbed (irregular and asymmetric) galaxy. We cannot detect the 
fainter arms very clearly in the Pa-a data. However, the central region surrounding the brightest 
two star-forming knots does appear kinematically disturbed - lacking structure in velocity and 
velocity dispersion. 

While a face-on disk could explain the absence in velocity structure, lack of structure in a 
makes that scenario less plausible. Therefore, inclination is unconstrained. Using the integrated 
spectrum, we find a =86 ±3 km s With no observed velocity gradient, we can only place 
an upper limit: v/cr <0.7, using Vc = (vmax— Vmin)/2 and no inclination correction. This case 
represents disturbed kinematics with multiple cores and underlying faint features, suggestive of 
recent merger activity. 



3.3. ScUVLG210358 

While morphologically similar to ScUVLG143417, we find that the kinematics in this galaxy 
are quite different. Again, the optical image reveals faint, irregular and asymmetric arms with 
a bright, unresolved star-forming core. However, in this case, there is a single bright, symmetric 
core. This galaxy is one of the most massive and luminous ScUVLGs in our sample. The position 
angle for the Pa-a isophotes (from the S/N contours) does not align with that of the velocity field, 
possibly indicating the presence of a warp or bar. The velocity gradient across the highest S/N 
region of this galaxy is very regular, maintaining similar structure to the outer, fainter regions of 
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the system. Hints of a merger are apparent from HST images, displaying irregular, clumpy light 
distribution and a potential faint companion to the south-west. Yet the dynamics of this galaxy's 
central region show ordered rotation, with the a peak corresponding with the galaxy's nucleus. 
Therefore this galaxy has dynamics consistent with disks. We find an inclination of i = 51 ± 4°, 
circular velocity of Vc = ^^^^ini =186^g^ km s~^, and a =175 ±4km s^'^, giving v/a ?al.l. Using 
Eqn. m at R=1.5kpc, we derive Mdyn = 4 x lO^^M© (for C=3.4) and exlO^^M© (for C=5). Similar 
to ScUVLG092600, the signature for rotation is present, yet, this system is marginally dispersion 
dominated. 



3.4. Simulation to high-z 

To best compare our results with high-z IFS observations (tracing Ha emission in the Kg- 
band), we perform artificial redshifting to determine how our 3 galaxies would appear if placed 
at z ^2. Using the recovered flux, velocity, and dispersion maps described in §3 we generate a 
synthetic model data cube which is rebinned to the angular scale which would be observed at redshift 
z = 2.2645 (chosen so that Ha emission in Ks band avoids strong OH emission features). Assuming 
that Ha emission traces Pa-a, we normalize this model to the total Ha luminosity determined 
by the unresolved SDSS s pectroscopic data and artifi cially 'observe' the model galaxy using the 



O SIRIS simulation code of lLaw. Steidel. &: Erbl (120061 ) and data reduction and analysis techniques 



of lLaw et al.l (|2007l ). We refer the reader to these papers for a full description of the routines, noting 



that the simulation code fully accounts for the characteristics of the OSIRIS+L GSAO system an d 



has been shown to well represent actual IFS observations of high-redshift galaxies (jLaw et al.ll2007l ). 
Each of our 3 galaxies is artificially redshifted and re-observed, using a simulated 50 mas lenslet 
scale in Ks band for a total of 2 hours of integration time under typical atmos pheric cond i tions . 



Thes e parameters are chosen to represent typical observational values used in iLaw et al.l (|2007l . 



20091 ). 



We show the high-z simulated data in Figure HI with black contours marking the z -^0.2 
data. For example, the middle panels display high-z velocity with low-z velocity (see data in 
Figure [2|) contours overlaid. In all cases, the high-z data matches the brightest low-z features. 
However, this figure illustrates some of the difficulties of high-redshift LEG observations: resolution 
loss in the simulated S/N data causes ScUVLG092600 to appear more face-on; the simulated 
(T-map of ScUVLG092600 appears more disturbed as edge features dominate; lower resolution 
blends the double core of ScUVLG143417 into a single blob; selecting velocities where S/N exceeds 
detectable limits gives the velocity map in ScUVLG143417 a slightly more regular appearance; 
surface brightness dimming hides the complex fainter structure in ScUVLG210358. We derive 
values for axis ratios, inclination, circular velocity, a and v/a from the high-z simulated data 
(shown in Table 1, alongside the observed low-z results). 

The v/a for high-z simulations of ScUVLG092600, ScUVLG143417 and ScUVLG210358 are 
l-Olgei 0.4zb0.1(no inclination correction), and 0.9±0.1, respectively. Therefore in all the cases, v/a 
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in high-z simulations differs somewhat from low-z values (as shown in Table 1). For ScUVLG092600, 
v/o" increases in the high-z simulated case, while decreasing for the other two galaxies. We find 
that the high-z simulations physically resen ible and have sim ilar, though slightly higher, v/o" values 
as z =2 LBGs, which have v/a ~ 0.3-0.7 (|Law et ahlbood ). 



4. Discussion and Summary 



We are pioneering a new study of the dynamics and kinematics in local (z <0.3) UV-selected 
galaxies that resemble LBGs in most physical properties. Similar studies have yielded mixed results 
about LEG dynamics at high-z, but our study of 3 low-z ScUVLGs shows that dispersion dominated 
kinematics is a viable scenario. Although 2 of the 3 cases show some signs of smooth rotation, 
v/fj < 1 in these galaxies. Since HST images show disturbed morphologies, nearby companions, 
and evidence of ongoing interactions (tidal features, plumes, etc.), it is not surprising that velocity 
dispersions in these galaxies are so high. 

IFS studies of dynamics in LBGs draw conc lusions ranging from smooth rotation to disper- 



Genzel et al.l (|2008l ) study a sample of massive (M 



10^ 



sion dominated kinematics, 
z = 2 — 3 rest frame UV-selected galaxies and classify these as turbulent, clumpy disks with formed 
bulges, undergoing secu l ar evo lution. Studying a subset of LBGs having the highest velocity fields, 
Forster Schreiber et al.l (|2006l ) find that the majorit y exhibit smooth rotation , consistent with cos- 
mological predictions of flattened disk- like galaxies (jMo. Mao. &: Whitdll998l ). while a more com- 
prehensi ve sample demonstrates larger c ontributions from both majo r mergers and high dispersion 
systems (IForster Schreiber et al.l l2009l ) . iNesvadba et al.l (j2006l ) and IStark et al.l (120081 ) determine 
gravit ationally lensed LBG s at z =3.2 and 3.07, respectively, are examples of rotationally supported 
disks. IWright et al.l (120071) fit their z=1.5 LBG with a rotation model but cannot exclude a merger 
scenario. iLaw et al. ( j2009l ) argue for largely dispersion do minated kinemati c s in t heir sample of 
LBGs. Analyzing IFS data for one massive LBG at z =3.2, INesvadba et al.l (|2008l ) find evidence 
for dispersion dominated kinematics. 

Comparing our sample of LBG analogs to LBGs at hi gh-z, our re s ults a re more consistent with 
compact, kinematically distur bed systems (i.e. similar to iLaw et al.l (120091 ). and see discussion b y 



Forster Schreiber et al 



Forster Schreiber et al 



(2009)) rather than larger, more massive disks [i.e. iGenzel et al.l (|2008l ). 
(120061 )). While our selection criteria for ScUVLGs includes compactness 
of the UV-emitting region, HST images do reveal underlying low surface brightness, extended 
structure in two cases ~ ScUVLG1434 and ScUVLG2103. These two examples represent the most 
massive ScUVLGs. Meanwhile, for high-z observations, optimization of candidates may bias this 
study by selecting larger and brigh ter galaxies and may i iot represent the g eneral physics governing 
galaxy behavior at these redshifts (jNesvadba et al.ll2008l : IStark et al.ll2008l ). Our view at low-z may 
provide insight about smaller, more compact, 2 = 2 — 3 LBGs. In these low-z cases, we argue that 
mergers may be responsible for the observed kinematics and morphologies and speculate that a 
similar scenario driving LBG characteristics at high-z cannot be excluded by current data (see also 
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Overzier et all (12003)). 



Lehnert et alj (|2009l ) investigate other mechanisms capable of causing high dispersions in high- 



z star-forming galaxies. They find that gas accretion cannot sustain the observed high dispersions. 
Studying the relationship between velocity dispersion and emission line characteristics, they con- 
clude intense star formation can account for high dispersions in these systems through a feedback 
cycle between turbulence in the interstellar medium (ISM) and star formation. The proposed sce- 
nario may offer a plausible expl a.nation for observed gas kinematics once galaxies reach a turbulent 
phase and associated starburst (jLehnert et al.l 12003 ) . 



We show that simulated high-z data produces v/o" values somewhat different from actual 
values. Misinterpretation of kinematic features may result from surf a,ce brightn e ss ari d resolution 
effects. However, these simulations do resemble z =2 LBGs studied by lLaw et al.l (|2009l ) in physical 
appearance and v/cj. 

Having launched this pilot study, we will cont inue researching a larger sample of ScUVLGs - 
performing kinemetry analysis (jShapiro et al.ll2008l ) in order to draw statis tically robust conclusio ns 
about the dynamics of these low-z, high SFR, UV-selected galaxies (see Goncalves et al. (j2009l )). 
Detailed IFS studies of local LEG analogs offer a unique method for capturing important physics - 
pinpointing the mechanisms most generally responsible for the high star-formation phase occurring 
in LBGs at high-z and in ScUVLGs at low-z. 



We thank the referee for helpful and insightful comments. The authors wish to recognize 
and acknowledge the very significant cultural role and reverence that the summit of Mauna Kea 
has always had within the indigenous Hawaiian community. We are most fortunate to have the 
opportunity to conduct observations from this mountain. A.R.B. wishes to thank Dr. Jacqueline 
van Gorkom for her advice. 
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Table 1. Summary of Observations 
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ScUVLG092600 141.502 44.460 0.181 0.101(0.208) 
ScUVLG143417 218.571 2.128 0.180 0.098(0.197) 
ScUVLG210358 315.995 -7.467 0.137 0.065(0.144) 


8 
15 
68 


9.20^^ 
10.89^= 
11.05^* 


0.78±0.05 
0.65±0.05 


40±5 
51±4 


31I3 84±3 
60 86 ±3 
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<0.70 
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0.79±0.02 


26±4 
39±2 


59t^^ 

35 
167^7 


59±17 
98±12 
175±11 


1 +"-^ 

^■^ -0.6 

<0.35 
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'^SpsF gives observed FWHM of the PSF, spatially smoothed FWHM is given in parentheses. 

''SFR was determined from 1.4GHz-derived SFR - proxy for total, unattenuated SFR (see lSasu-Zvch et al.l |20o3)) 
'^M* based on SED-fitting Csee lSalim et al.l l|2005l ')'). 

'^M* is dust corrected median stellar mass llKauffmann et al]|2003f) from MPA/JHU value added catalog (see |http: //www.mpa-garching.mpg.de/SDSS/index_old.html^ . 
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Fig. 1. — HST optical (red and green) + UV (blue) images (left to right, see lOverzier et alJ ([200a, 
2o3)): ScUVLG092600, ScUVLG143417, and ScUVLG210358. Magenta boxes mark OSIRIS FOV. 
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Fig. 2.— OSIRIS IFS results for three ScUVLGs- ScUVLG092600 (top), ScUVLG143417 (middle 
row) and ScUVLG210358(bottom). We show S/N maps (left), velocity maps (middle column), 
velocity dispersions (right) with S/N contours overlaid and nucleus (center of the brightest region) 
displayed as a white triangle. Contour levels correspond with marked values of the S/N color bar. 
The smoothed PSF and physical scale are shown at upper left. 
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Fig. 3. — PV diagram for ScUVLG092600 is shown at bottom for a slice along the velocity gradient 
(marked in magenta, on the velocity field, at top). PV diagram illustrates slight velocity gradient 
with high dispersion. 
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Fig. 4. — High-z simulations for S/N, velocity, and cr maps (from left to right column) for ScU- 
VLG092600 (top), ScUVLG143417 (middle row), and ScUVLG210358 (bottom). Overlying black 
contours outline data from Fig. [21 such that contours show low-z velocity data in middle column 
and low-z a data in the right column, with levels matching marked values from color bars in Fig. 
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